Exponential Convergence of Cellular Dynamical Mean Field Theory: Reply to the 
comment by K. Aryanpour, Th. Maier and M. Jarrell (cond-mat/03014603) 
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We reply to the comment by K. Aryanpour, Th. Maier and M. Jarrell (cond-mat/0301460 ) on 
our paper (Phys. Rev. B 65 155112 (2002)). We demonstrate using general arguments and explicit 
examples that whenever the correlation length is finite, local observables converge exponentially fast 
in the cluster size, L c , within Cellular Dynamical Mean Field Theory (CDMFT). This is a faster 
rate of convergence than the behavior of the Dynamical Cluster approximation (DCA) thus 

refuting the central assertion of their comment. 



The development of cluster extensions of dynamical 
mean field theory is an active area of research. Cluster 
dynamical mean field theories, construct approximations 
to the solution of model Hamiltonians on the lattice in 
terms of the solution of a cluster impurity model. Differ- 
ent cluster schemes should be viewed as different trun- 
cations of the full quantum many body problem. As the 
size of the cluster tends to infinite, all cluster schemes 
approach the exact solution of the lattice model. A rele- 
vant question is, for a given computational power, (which 
only allows the investigation of small cluster sizes) which 
truncation is closer to the result in the thermodynamic 
limit. 

Our previous publication! investigated and extended 
two cluster schemes, the cellular dynamical mean 
field theory^ (CDMFT), and the dynamical cluster 
approximation 3 (DCA) by applying it to an exact solv- 
able model and concluded that CDMFT converges faster 
than DCA to the exact solution of that model. In their 
comment 4 , Aryanpour et. al introduce a new general- 
ization of the original DCA equations that take into ac- 
count better the non-local interaction, and they argue 
that their new method converges faster than CDMFT to 
the exact solution of the model as the size of the clus- 
ter increases. They also comment that our findings of 
rapid convergence of CMDFT are surprising, in the light 
of an earlier publication^ in which they concluded that 
CDMFT converges to the infinite cluster size limit with 
corrections of order 0(1/ L c ) where L c is the size of the 
cluster while DCA converges faster, with corrections of 
order 0(1/1%). 

In this reply to their comment we point out that local 
observables in CDMFT generally converge exponentially 
at finite temperatures, as long as the relevant correlation 
length is finite (a situation that excludes a critical point). 
This statement persists at zero temperature in systems 
which have an energy gap. 

We demonstrate the exponential convergence of CDMFT 
in three steps. First we present general arguments in 
favor of exponential convergence of local observables in 
CDMFT whenever the relevant correlation length is fi- 
nite. This is a direct consequence of the cavity construc- 



tion underlying the method. In the process we explain 
why the convergence criterion introduced in reS and used 
to conclude that CDMFT converges as 1/L C is not an 
appropriate measure of convergence of local observables 
in CDMFT, which instead converge much faster than a 
power law in L c . Then we present the numerical results 
for the SU(N) spin chain studied ini that agree com- 
pletely with our general arguments and we display ex- 
plicitly an example of the exponential convergence in this 
model. Finally, we discuss another simple case, the semi- 
classical limit of the Falikov-Kimball model in one dimen- 
sion. Previous work established 7 that in this limit quan- 
tum cluster approximations reduced to classical cluster 
approximations. This allow us to compare DCA and 
CDMFT in detail using simple analytical considerations. 
In particular we unveil, in an explicit example, that the 
DCA predictions for local observables converge as 1/L 2 C 
even when the same quantities obtained solving a finite 
system of size L c (with for example periodic boundary 
condition) converge exponentially fast (in L c ) to their 
thermodynamic limit value. 

Let us start with some general considerations. If one 
was able to trace out exactly all the degrees of freedom 
outside the cluster to get an exact effective action for 
the degrees of freedom inside the cluster then the trans- 
lation invariance of the effective action would be broken 
(degrees of freedom near the boundary of the cluster are 
affected by the enviroment more than the bulk degrees 
of freedom) but observables within the cluster would still 
be translationally invariant. So on very general grounds 
we expect that the Weiss field, (or hybridization function 
A) which describes the effects of the degrees of freedom 
integrated out, is large near the boundary and small (in 
fact exponentially small if the correlation length is finite) 
inside the cluster. 

CDMFT is an approximate way to realize this cavity con- 
struction. It produces a Weiss field which is large at the 
boundary and small in the bulk, but because of its ap- 
proximate character it produces non translation invariant 
observables. However, whenever the correlation length 
is finite: (1) bulk quantities for a free system of size L 
with, for example, free boundary conditions converge ex- 
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ponentially fast in L to their thermodynamic limit (2) the 
CDMFT approximation should improve the convergence 
of bulk quantities, which as a consequence should be at 
least as fast as the one of the system with free bound- 
ary conditions. Note that this is not the case of DCA 
that still converges as 1/L% worse than the results for a 
finite system with, for example, free boundary conditions 
which would converge exponentially fast in L c \ We will 
discuss an explicit example of this behavior below. 

Because CDMFT breaks the translation invariance in- 
side the cluster it is important to extract properly the 
value of local observables. These are well represented in 
the center of the cluster (see Fig. 1) , and will have more 
deviations near the boundary. The role of the cavity field 
is to try to reduce these deviations, but in approximate 
treatments there will always be errors of order one near 
the boundary. Aryanpour et al^ concluded the 1 / L c con- 
vergence by estimating the value of a local observable of 
interest by doing a flat average (of the value of the local 
observable) over the cluster. Since the error is of the or- 
der one on the boundary, one obtains an error that dies 
out as the ratio surface over volume of the cluster, i.e. 
1/L C . On the other hand, if we extract this value from 
the center of the cluster, which is the natural prescrip- 
tion dictated by the approximate CDMFT cavity con- 
struction, we obtain exponential convergence as we show 
below. 

We now turn to the simple one dimensional SU (N) chain 
studied in 1 whose Hamiltonian is: 

H = (!) 

2iV ^ ' (ft&f i ' a ''fi+l,cr'f i + 1 - a ' fi+l,crf i + 1 ,< J 'fi,a'f i >< J ) 

where / are fermions operators, i = 1, . . . , L and a = 
1,...,N and we take the large L and N limits. This 
model is a generalization introduced by I. Affleck and B. 
Marston& of the Hubbard-Heisenberg model where the 
SU(2) spins are replaced by SU(N) spins, the on site re- 
pulsion is scaled as 1 /N and the large N limit is taken. 
One can apply CDMFT and DCA to this model. Note 
however, that because the interaction is non-local, there 
are different possible extensions of usual cluster methods 
to this case. We extended DCA in a way based on the 
real space perspectivei, Aryanpour et. al introduced a 
different generalization of DCA which takes into account 
better the non-local interactions. Our procedure is there- 
fore not an incorrect application of DCA, as claimed by 
Aryanpour et al., but only a different generalization of 
DCA to the case of non-local interaction. The results of 
the two different generalization are discussed in 1 *^ and 
in the following we shall focus on the generalization of 
Aryanpour et al. which has been shown to converge to 
the thermodynamic limit with an error of the order of 
1/L% where L c is the size of the cluster. This rate of 
convergence is a general property of DCA^, at least far 
from critical points. 



In the following we will use J as the unit of temperature 
and therefore we put J — 1 and we rescale the hopping 
term t — > t/J. The thermodynamics of this model can 
be solved exactly since in the large N limit the quantity 

X = jj J2a /i,o-(*)/i+i,o-(*) does not fluctuate. Indeed (jTJ) 
reduces to a free-fermions Hamiltonian with a "renor- 
malized" hopping term t —> t + x an d a self-consistent 
condition on x' 

X = ~ f^ E k) cos fc > E k = -2(t + X) cos k + n 

k 

(2) 

where /i is the chemical potential, f((3Ek) is the Fermi 
function and (3 is the inverse temperature 1 . 
DCA and CDMFT result in self-consistent eqs. for 
Xx = jjYscr < fx,afx+i,o > 1A - Because of the trans- 
lation invariance of DCA the Xx CA are independent of 
x inside the cluster. Once self-consistency is achieved 
one can use the cluster quantity to obtain the DCA 
lattice prediction Xum A as explained in 4 . This quantity 
converges to the exact x with corrections of 0(1/ L^)i£. 
Contrary to DCA, the self-consistent CDMFT eqs. 
breaks the translation invariance inside the cluster. 
In fig. n we plot the error of the CDMFT prediction on 
X: Xx — x, as a function of x, where Xx is the CDMFT 
solution for /3, the inverse temperature, equal to 4, 
t = fi = 1 and L c = 38. 
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FIG. 1: Xx — X as a function of x, for 0, the inverse temper- 
ature, equal to 4, L c — 38 and t = fi = 1. 

This figure clearly shows the behavior discussed previ- 
ously, namely translation invariance inside the cluster is 
broken: errors are smaller in the bulk while they remain 
of the order one at the boundary. Aryanpour et al^ con- 
cluded, by carrying out a flat average Xx/ {L c —1) over 
the cluster, that the error within CDMFT is expected to 
be of the order 1/ L c (more generically is surface over vol- 
ume, hence, 1/L C also in dimension larger than one). As 
discussed above, it is better to extract the CDMFT esti- 
mators weighting bulk values more than boundary values. 
The easiest thing to do is just taking the value of Xx a t 
the center of the cluster. In fig. [21 we compare the er- 
ror obtained doing the flat average (square, dotted line) 
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to the one obtained focusing on the bulk values (circle, 
solid line). This figure conveys two important informa- 
tion: first the Axtuik = \xl c /2 ~x\ 1S much smaller than 
the flat average one Ax/ a = \J^ x Xx/(L c — 1) — x| for 
large values of L c . Second, as shown in the inset, the 
error Axbuik multiplied by L\ is still decreasing fast as 
a function of L c , i.e. Axbuik decreases much faster than 
1/i 2 . . Instead the error corresponding to the flat av- 
erage leads to a straight line corresponding to the 1/L C 
behavior discussed above and irA Note that in this plot 
the DCA prediction would lead to a function approach- 
ing a constant when L c — > oo. Finally, in fig. [21 we plot 




FIG. 2: Axbuik (circles, solid line) and Ax/ a (squares, dotted 
line) as a function of L c . Inset: AxbuikL^ (circles, solid line) 
and AxfaLc (squares, dotted line) as a function of L c . 

A-Xbuik m a logarithmic scale as a function of L c . The 
exponential convergence is manifest (the straight line is 
a guide for the eye) and in complete agreement with our 
general discussion above. 



FIG. 3: Axbuik as a function of L c in a logarithmic scale. 

Another simple example that sheds light on the con- 
vergence properties of CDMFT and DCA is the Falikov- 
Kimball model in the large U limit. In this limit the 
system becomes a classical Ising model and the quantum 
cluster schemes map on classical cluster schemes^. In the 
following we will focus on the easy one dimensional case 



in the paramagnetic phase. In this case, since the Weiss 
field is zero, doing CDMFT on a cluster of size L c cor- 
responds only to solving a finite system of size L c with 
free boundary conditions. Instead DCA corresponds to 
solving a finite system of size L c with periodic boundary 
conditions and a coupling 
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where J is the original spin-spin coupling. So, DCA 
and CDMFT differ in the different boundary conditions 
but also in the fact that for DCA the internal spin-spin 
coupling is modified from its original value whereas for 
CDMFT is not. This leads to very different convergence 
properties. Let us focus for example on the prediction for 
the spin-spin correlation C =< SiS i+ i >. The CDMFT 
and DCA results read: 

Ccdmft = tanh/?J 



C 



DCA 



t&nh(3J D cA + (tanh/? Jdca) 1 



1 



(tanh f3 J ' DC a) l " 



In this case the CDMFT prediction is exact because trac- 
ing out the spins outside the cluster indeed leads to a zero 
magnetic field on the boundary but this is of course a pe- 
culiarity of this simple case. Instead, there are two types 
of corrections to DCA (note that Jdca = J +0{1/L C 2)). 
There are corrections which die out exponentially fast 
as oc (tanh/?J) Lc . These are the same type of correc- 
tions that one obtains using periodic boundary conditions 
or other type of boundary conditions for a free system. 
However, there is a much larger correction coming from 
the first term in the numerator of Cdc a which leads to 



Cdca = tanh /3 J — 



/?Jtt 2 



3(cosh/?J) 2 (2L c ) 2 



+ 0(L-*) 



Thus applying DCA to this problem one obtain a 
convergence as 1/L 2 . which is much worse than the 
exponential one corresponding to solving the free model 
with periodic boundary conditions. The origin of this 
behavior can be traced to the fact that as DCA forces 
translation invariance inside the cluster the couplings 
are changed everywhere in the system by an amount of 
the order 1/1%. So even if the correlation length is finite 
this error dominates the convergence. 
While we stress the obvious advantages of CDMFT, it is 
also worthwhile to point out the aspects of the CDMFT 
method (and cluster methods in general) that still 
require development. The lack of translation invariance 
of CDMFT , which in the toy model manifests itself 
in the site dependence of the bond expectation value 
Xx is certainly one of them. For example, CDMFT 
predicts a finite temperature phase transition for the one 
dimensional Falikov-Kimball model in the semiclassical 
limi1i£i£. This is due to the fact that value of the 
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Weiss field on the boundary is strongly coupled to 
the value of the propagator on the other boundary 
and, unfortunately, as discussed previously, the error is 
much larger at the boundary than in the bulk. This 
may not be a serious problem for phases with broken 
symmetry, but certainly is in translationally invariant 
phases. One possible solution of this problem would 
be to modify the self-consistent equations that express 
the Weiss field as a function of the propagator so as to 
use more heavily bulk values of the propagator, which 
have a reduced error relatively to the boundary. Clearly, 
further investigations are needed to optimize CDMFT 
in the light of this point. A related issue, stressed in 
refii^, is that the lattice self energy in CDMFT is a 
derived quantity, obtained from the cluster self energy 
entering the CDMFT equations. The lattice self-energy 
is obtained using a matrix w, via the formula (10) in 
re£k If w is positive definite, then the lattice self-energy 
is causal, but one can sometimes obtain better estimates 
by using other matrices. For the toy model, the matrix 
w are not restricted since the self energy is real. The 
positive definitiveness of w is a sufficient, but not a 
necessary condition to maintain causality. A general 
constructive way to find the best estimator of the lattice 
self energy (preserving causality) for an arbitrary model 
Hamiltonian is lacking. However, an important criterion 
to follow is trying to extract the prediction on the lattice 
self energy from the bulk values of the cluster self-energy 
in order to minimize the error as discussed previously. 
None of these issues, however, affect the basic fact that 
expectation values of physical observables which are 
local (namely defined on a restricted neighborhood in 



physical space) converge as l/L% for DCA as shown 
ini, while for CMDFT they converge exponentially 
in situations like the one outlined for the toy model, 
when the correlation length is finite, and not as 1/L C 
as claimed by K. Aryanpour et ali By exploiting 
the freedom in the choice of basis, which is inherent 
to the original CMDFT formulation, one can improve 
convergence of observables which become local when the 
approximation is formulated in a different basis set. The 
problem of convergence as function of cluster size at 
zero temperature or at a quantum critical point, or for 
quantities that are dominated by massless excitations, 
remains an open question. However these problems 
could be better addressed by techniques other than 
quantum cluster methods. 

Finally, most studies can only be done for small clusters, 
and it is important to understand whether the results 
obtained in small clusters are representative of the 
thermodynamical limit. Recent CDMFT studies of the 
Hubbard model, in one dimension, show that while 
even - odd effects are important, even clusters of small 
size can give very accurate results^ as compared with 
exact Bethe Ansatz results in the thermodynamical limit. 
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